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Abstract 

A simplified method for stratified floor strata ultimate bearing capacity (UBC) 
assessing, in room-and-pillar mining conditions, is presented. The physical problem 
consists of the following two different cases of rectangular rigid pillar of size BxL 
resting on a two-layer floor system: (a) with the upper weaker layer of thickness H 
overlying a stronger, infinite rock deposit. It has been assumed that bearing capacity 
of floor strata may be analyzed as a shallow foundation problem, using the general 
bearing capacity equation given by Brinch Hansen with appropriate shape, inclination, 
and surcharge depth factors, or may be treated as a punching failure problem, 
particularly in a case of relatively strong upper stratum presence. A broad review of 
available analytical techniques did show that a finite thickness of the upper weaker 
layer causes a significant increment/decrement in its ultimate bearing capacity in 
comparison with homogeneous conditions. The proposed approach has been also 
validated by numerical modeling utilizing FLAC3D the finite difference computer code. 

Introduction 

Presently, design of hard rock pillars for weak floor strata conditions is 
based on the weak floor strata ultimate bearing capacity (UBC). Design 
techniques utilized here typically estimate UBC underneath full size pillars 
using Vesic [1] analysis which considers weak floor strata as a two-layer 
cohesive soil or rock system, with a weak layer overlying a stiffer layer and 
the angle of internal friction {§) for both layers equaling to zero. The 
cohesion (c) of the weak layer is estimated more often from its moisture 
content [2] than from the laboratory compressive strength data as proposed 
by Vesic. This technique, called Vesic-Speck approach, is used most 
commonly for mining applications in the U.S.A. today; however, this analysis 
technique has never been validated for shallow foundations on layered rock 
strata involving a weak layer(s) overlying a stronger rock layer. To overcome 
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most of those limitations, the analytical Pytel-Chugh's approach has been 
developed [3], including the effect of adjacent pillars and the non-zero values 
of the angle of internal friction for both layers. Since then this technique has 
been successfully applied mainly in Midwestern coal mines in U.S.A. 

It has been assumed that bearing capacity of coal pillars on floor strata 
may be analyzed as a shallow foundation problem, assuming the validity of 
the foundation general bearing capacity. All field observations show that a 
finite thickness of the upper weaker layer causes a significant difference in 
its ultimate bearing capacity in comparison with homogeneous conditions. 
For a purely cohesive floor system {§=0), the Vesic's approach is most 
usefull and convenient; but where cohesion and the angle of internal friction 
for weak floor strata can not be assumed to be equal to zero, the Mandel 
and Salencon [4] single-layer technique, based on the slip-line method with 
Coulomb criterion of failure is more suitable. 

Figure 1 depicts the problem where multiple pillars of width B and length 
L, spaced s apart, are resting on a finite layer of of weak floor stratum (ci, fa, 
Yi) of thickness H underlain by a deformable infinite layer with parameters c 2 , 

fa, Y2- 
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Fig. 1 . Physical Problem 


It should be noted that currently updated Pytel-Chugh's method is the 
only technique which considers the effect of adjacent pillars and the angle of 
internal friction not equal to zero for both floor deformable layers. To date, 
the appropriateness of this design technique has been also established 
through field observations in Polish deep copper mines [5]. This approach 
may be summarized as follows: 

- the foundation is considered as consisting of two layers, with a weak 
floor layer overlying the stiff layer, 
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the effect of adjacent pillars and non-zero values of the angle of 

internal friction for both layers is included, 

a new modified bearing capacity factor Nt is calculated based on a 

combination of two theoretical solutions: the Mandel-Salencon [4] 

bearing capacity estimation approach (effect of finite thickness of the 

upper stratum), and the Mandel [6] approach (effect of the pilars' 

spacing). 


1. The Proposed Approach 

The ultimate bearing capacity of two-layer floor system q (Fig. 2) may be 
expressed by the following basic approximation: 


- = — + — ,or q= W1W2 , [q,<q 2 ) 
q q, q 2 q^+q 2 
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Fig. 2. Proposed approach for calculation of pillar bearing capacity 
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where: q 1 = q 1 -=^- represents the UBC of pillar substrata, assuming the lower 

layer to be infinitely rigid {F cV F c2 Mandel & Salencon coefficients, Fig. 3) and 
spacing s between pillars is greater than: 


s = 2 -"— cos(^ - %) exp[(^ - 1) tan * ] (2) 

cos^ 4 2 4 2 

q 2 =^q 2 represents bearing capacity of the lower layer including factor of 

s 
majoration: £(— ,<f> 2 ) due to the presence of adjacent pillars [6], (see Fig. 4), 
B 

while q v q 2 axe the ultimate bearing capacities of the floor strata composed of 

layers 1 or 2 exclusively (homogeneous conditions), subjected to the load 
transmitted from the pillar: 

9i = CiA/ c1 s c F c1 + - yfiN^s^ , q 2 = c 2 N c2 s c F c2 + - r 2 BN y2 s r F r2 

where: F c and F r are the Mandel & Salencon's coefficients for homogeneous 

conditions for a given value of angle of internal friction tyn (for § = 0° - F c = 

0.778, for <|> = 10° - F c = 0.691 ; for § = 20°, F c = 0.585; for ty = 30°, F c = 

0.460), N c and N y are bearing capacity coefficients and s c are shape factors. 

Actually there is no real physical basis for Equation 1 . It is a kind of formal 
tool which couples two separate analytical solutions since there is no closed 
expression for the general two-layer floor problem. The proposed form of 
Eq. 1 is the simplest one which satisfies exactly the mentioned boundary 
conditions, however for finite values of qi and q 2 it produces an error which 
can not be determined a priori. Therefore this equation has been 
transformed into the following form: 

_q 1 q^ + E %_ t ( < ) (3) 

9i + 9 2 9 2 

where E r is the maximum error which occurs when q^ =q 2 . From Eq. 3 one 
may obtain: 


F =92. 
9i 


Q _J7i92_ 
9i+9 2 


(4) 
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The error term is assumed to have the form E. ^because it is able to 

express conveniently the deviation from the exact value in function of 
difference between two layers' strength parameters. If q., = q 2 then Ci = c 2 

and §i = §2 and then two-layer system transforms into one-layer system 
(independently on H/B) with row of pillars on the boundary surface and then 
Eq. 4 may be simplified into the following form: 


<7i 


<f- 


1 


1 1 

— + — 


(5) 


where: F 


-, F c is the Mandel-Salencon coefficient (Mandel & Salencon, 


1969) suitable for angle of friction between coal pillar and the upper floor 

2 
stratum S = —fa (values of F c may be find in Figure 3), £is the Mandel factor 

of majoration (Mandel, 1965) calculated using a given distance s and angle (j) 2 
(values of £are presented in Figure 4). Substituting (5) to (1) and assuming a 
weightless material underneath the shallow pillar we finally obtain: 


q=q,N t = q, 
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(6) 


2. Numerical Example 

Design approaches presented in the paper are illustrated with a physical 
problem solution with the following geometry: 


- ratio of pillar length over its width: L/B = 1 00, 

- ratio of opening width over pillar width: s/B = 4.0 (for results shown in 
Fig. 5) and s/B = 0.5 (Fig. 6) 

- angle of internal friction for floor strata: (j) 2 = 40°. 
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Fig. 3. Values of F c of Mandel-Salencon (1969) Coefficient F c 
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This physical problem was also analyzed using different values for relative 
floor strata thickness H/B and the ratio -^-. Calculation results are 

9 2 

presented in Figures 5-6. 

Differences between results obtained using the proposed approximate 
approach and the finite difference method (FDM) are presented in Fig. 7 for 
a chosen case of underground geometry and geotechnical conditions. 

Conclusions 

1. The updated analytical expression for UBC of immediate two-layer weak 
floor strata (Eq. 6) is more realistic than its previous version since it has 
been associated in some instances (e.g. for very low H/B and high values 

of § 2 or for -=- close to 1.0) with the unreasonable estimates of q greater 

9 2 

than q 2 . 

2. The effect of pillar spacing is very important particularly for lower values of 
H/B. 

3. The effect of ^- ratio on the resultant value of ultimate bearing capacity q 

Q 2 
is almost linear in a case where the lower, stronger stratum is located at 

the moderate depth (H/B =0.2) -> q = ^-q 2 . 

9 2 

4. The updated estimate of immediate floor strata bearing capacity has been 
proved to be a very convenient and relatively accurate analytical tool for 
room-and-yield pillar systems design. 
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Fig. 5. The effect of H/B on two-layer immediate floor strata UBC 
d = 0.5 MPa (- - -), ci = 2.0 MPa (— ), c 2 = 2.0 MPa, 4> 2 = 40° , s/B = 4.0 
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Fig. 6. The effect of 3- ratio on UBC of immediate floor strata 
Ci = c 2 .= 2.0 MPa, fa = 40° , s/B = 0.5 
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Wytrzymatosc spqgu uwarstwionego obciazonego pod filarem 

Slowa kluczowe: staba warstwa spqgowa, opor graniczny podloza, metoda roznic 
skohczonych 

W artykule przedstawiono uproszczonq metode szacowania nosnosci 
uwarstwionego spqgu, dla warunkow eksploatacji metodq filarowo-komorowa. Model 
fizyczny zagadnienia opisuje przypadek prostokqtnego, sztywnego filara 
o wymiarach BxL spoczywajacego na dwuwarstwowym ukfadzie spqgu 
bezposredniego, ze slabszq warstwq potozonq wyzej o grubosci H, zalegajqcq nad 
mocniejszq warstwq skalnq o nieskohczonej miqzszosci. Przyjeto, ze nosnosc 
podloza mote bye analizowana jako przypadek fundamentu bezposredniego 
poprzez ogolne rownanie nosnosci zaproponowane przez Brincha Hansena, 
z zachowaniem odpowiedniego ksztattu, nachylenia oraz wspolczynnika 
zagtebienia. Przeglqd dostepnych matematycznych analiz stosowanych dla 
opisanego problemu wskazuje, ze skohczona miqzszosc slabszej warstwy powoduje 
wzrost/spadek wartosci nosnosci, w porownaniu z warunkami panujqcymi 
w gorotworze jednorodnym. Proponowane podejscie rozwiqzania problemu zostalo 
zweryfikowane przy zastosowaniu metod numerycznych wykorzystujqeych metode 
roznic skohczonych (F LAC 3D). 


